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ABSTRACT 

The aim of this study is to search for observational evidence of vertical iron stratification in the 
atmosphere of fourteen blue horizontal-branch (BHB) stars. We have found from our numer- 
ical simulations that five BHB stars: B22, B 186 in the globular cluster NGC 288, WF2-820, 
WF2-2692 in Ml 3 and B203 in M15 show clear signatures of the vertical stratification of iron 
whose abundance increases toward the lower atmosphere. Two other BHB stars (B334 in M15 
and B176 in M92) also show possible iron stratification in their atmosphere. A dependence 
of the slope of iron stratification on the effective temperature was also discovered. It is found 
that the vertical stratification of iron is strongest in BHB stars with T^g around 1 1,500K. The 
slope of iron abundance decreases as Teff increases and becomes negligible for the BHB stars 
with Teff ~ 14,000K. These results support the hypothesis regarding the efficiency of atomic 
diffusion in the stellar atmospheres of BHB stars with Tefl- ^ 11, 500 K. 
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1 INTRODUCTION 

According to the current understanding of stellar evolution, the 
horizontal-branch (HB) stars are post-main sequence stars that 
bum helium in their core and hydrogen in a shell (e.g. Hoyle& 
Schwarzschild 1955; Moehler 2004). In this paper we consider the 
HB stars that are located in the blue part of the HB, to the left of the 
RR Lyrae instability strip. These stars are commonly called blue 
horizontal-branch (hereafter BHB) stars to distinguish them from 
the red horizontal-branch (RHB) stars, which exhibit different ob- 
servational properties. 

BHB stars with Teff >: 11,500K show several observational 
anomalies. By studying the position of BHB stars of globular clus- 
ter^^on Hertzspruing-Russell diagrams, it is found that they exhibit 
photometric jumps (Grundhal et al. 1999) and gaps (Ferraro et al. 
1998) at Tgff ~1 1,500 K. These stars also show large abundance 
anomalies (Glaspey et al. 1989; Behr et al. 1999; Moehler et al. 
1999; Behr 2003a). For instance, iron is generally overabundant 
in these stars relative to their cluster abundance. It is also found 
that the spectroscopic gravities of these stars are lower than those 
predicted by the canonical models (Crocker et al. 1988; Moehler 
et al. 1995). It is commonly believed that atomic diffusion of the 
elements in the atmosphere of these stars is responsible for the 
aforementioned observational anomalies. Atomic diffusion arises 
from the competition between radiative acceleration and gravita- 
tional settling. This can produce a net acceleration on atoms and 
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ions, which results in their diffusion ( jMichaud 19701 1 and which 
may lead to abundance anomalies. 

In addition to the detection of vertically stratified abundances 
in some BHB stars (Khalack et al. 2007 and 2008, more informa- 
tion is given below), BHB stars with T^ff y 11,500K have lower 
rotational velocities than their cooler counterparts (Peterson et al. 
1995; Behr et al. 2000a and 2000b; Behr 2003b). These observa- 
tions suggest that the atmospheres of such stars are stable enough 
to have atomic diffusion. Theoretical results of Quievy et al. (2009) 
show that helium sinks in stars with low rotational velocities, which 
leads to the disappearance of the superficial He convection zone. 
This then opens the door for atomic diffusion to play a role. 

Other theoretical results also concur with the belief that dif- 
fusion is present in the atmospheres of BHB stars with T^ff >z 
11,500K. For instance, the atmospheric models of Hui-Bon-Hoa, 
LeBlanc & Hauschildt (20001 showed that the observed photomet- 
ric jumps and gaps for hot BHB stars can be explained by ele- 
mental diffusion in their atmosphere. These models calculate self- 
consistently the structure of the atmosphere while taking into ac- 
count the stratification predicted by diffusion (assuming equilib- 
rium). These models have recently been improved (LeBlanc et al. 
2009) and confirm that vertical stratification of the elements can 
strongly modify the structure of the atmospheres of BHB stars. 
Such structural changes for the atmosphere lead to the photometric 
anomalies discussed above. 

As mentioned above, by synthesizing spectral line profiles, 
Khalack et al. ( 2007 ) found vertical abundance stratification of sul- 
fur in the atmosphere of the field BHB star HD 135485. Vertical 
stratification of iron abundance was also detected by Khalack et al. 
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Table 1. Journal of Keck+HIRES 


spectroscopic observations of the selected hot BHB stars from Behr (2003a). 






Cluster/Star 


HJD 


Exposure 


S/N 


Coverage 


Tiff 


logg 


V sin i 


Vr 


[Fe/H] 


[He/H] 




2450000+ 


Time (s) 




(A) 


(K) 


(dex) 


(kms-') 


(kms-') 






M13A¥F2-820 


1052.7815 


3x1200 


34 


3888-5356 


11840 


3.90 


4.17 


-239.7±0.7 


-U.//_o.23 


I ng+0.47 
^■"^-0.44 


M13A¥F2-2692 


1047.7400 


3x1800 


34 


3885-6292 


12530 


4.08 


4.64 


-236.2±0.7 


-0.23+»]« 


1 01+0.96 


M13AVF2-3123 


1248.1298 


3x1200 


35 


3885-6292 


13670 


4.31 


6.32 


-238.3±1.1 


n 7, +0.65 


1 98+'-5' 


M13AVF3-548 


1053.7293 


3x1200 


32 


3888-5356 


13100 


4.16 


4.29 


-235.3±0.7 


-OMtlf, 


-2.io^i:§^ 


M13A¥F3-1718 


1248.0811 


3x1200 


48 


3888-5884 


11510 


3.79 


1.89 


-244.0±0.3 


+0.02;«;1« 


-i-86^^:S 


M15/B203 


1408.0180 


3x1200 


28 


3885-6292 


13990 


3.84 


4.89 


-94.7±1.2 


+0.02;«;1 


-2-54^i:^? 


M15/B315 


1409.0238 


3x1200 


24 


3885-6292 


12890 


3.81 


1.72 


-103.3±0.8 


+0.04;«;^ 


, 7^+0.94 


M15/B334 


1052.8791 


4x1500 


39 


3888-5356 


10750 


3.61 


9.22 


-109.3±1.3 


.9^7+0.19 
^■-''-0.21 


-U.D4 Q4J 


M15/B374 


1406.0183 


3x1200 


21 


3885-6292 


12820 


3.82 


3.92 


-106.5±0.9 


+U.24_Q37 


1 Q7+0.94 


M92/B176 


1049.0646 


3x1800 


56 


4020-5520 


11150 


3.76 


6.96 


-116.3±1.4 


-2.25+«i5 


n s5+°-3« 

-U.33_o32 


NGC288/B16 


1406.0850 


3x1200 


21 


3885-6292 


14030 


4.15 


3.08 


-44.7±0.9 


+0.66l«;« 


9 O9+0.91 
-2-02_o.71 


NGC288/B22 


1407.0955 


3x1200 


17 


3885-6292 


12130 


4.01 


2.48 


-43.8±0.9 


+0.28^U:i 


9 ^£+1.26 
^■^°-0.64 


NGC288/B186 


1409.0862 


3x1200 


31 


3885-6292 


11390 


3.94 


3.89 


-43.8±0.5 


+0.17^°:^] 


1 zl<;+0.65 


NGC288/B302 


1408.0891 


3x1200 


26 


3885-6292 


13230 


4.15 


1.85 


-43.8±0.5 


+o.6i^«:3^ 


9 ,-,+0.88 



( 12008 1) in the atmosphere of two BHB stars in the globular clus- 
ter M15 and for one star in the globular cluster M13. Also, more 
recently, Hubrig et al. ( 120091 ) have found isotopic anomalies of cal- 
cium in six BHB stars. 

In this paper we will attempt to detect signatures of vertical 
abundance stratification of iron in a set of fourteen BHB stars from 
line profile analysis for which we have appropriate data. This ex- 
tension of the studies of Khalack et al. (2007 and 2008) will serve 
to provide a larger amount of observational data for iron stratifica- 
tion which can then be used to verify theoretical models. We have 
selected the BHB stars from the spectra of Behr (2003a) relying 
on their value of effective temperature and their low Vsin/ (^ 10 
km s^') so that atomic diffusion phenomena would most probably 
be present and detectable in these stars. The properties of the ac- 
quired spectra are discussed in Sec.|2l while in Sec. [3] the details 
concerning the simulation routine and adopted atmospheric param- 
eters for the program stars are described. The evidence for vertical 
stratification of iron is given in Sec.|4]along with the estimation of 
the mean iron abundance, V sin / and radial velocities for each BHB 
star considered here. A discussion follows in Sec.|5] 



2 OBSERVATIONS 

The selected BHB stars were observed with the Keck 1 telescope 
and the HIRES spectrograph. The stars were selected from the large 
sample of stars observed by Behr (!2003ai The criteria by which 
the stars considered here were chosen is that their Teff is close to 
or larger than 11,500 K and that they possess a low Vsin/ (^10 
km s^'). These criteria were chosen to maximize the chance that 
atomic diffusion is present in their atmosphere and that its effects 
are potentially observable. 

Table [T] summarizes the information concerning the observa- 
tions of the fourteen stars chosen by giving (in individual columns) 
the object identification, the heliocentric Julian Date of the obser- 
vation, the exposure time, the S/N per pixel the spectral coverage, 
the size of the seeing disk (FWHM), the effective temperature, the 
logarithm gravity and the heliocentric radial velocity. 

The size of the CI slit is 0.86 arcsec which provides a spec- 
tral resolution ofR = X/&X =45,000. Behr ( |2003al ) has found that 
for the aforementioned stars the underfilling of the slit should not 
change the estimated spectral resolution by more than 4%-7%. 



Table 2. Properties of the BHB stars studied here. 



Cluster/Star 



[Fe/H] 



V sin i 
(kms-') 



Vy 

(km s" 



M13/WF2-820 


-0.09±0.32 


4.2±1.0 


-239.6±0.8 


79 


M13/WF2-2692 


+0.03±0.25 


4.4±1.0 


-236.2±1.0 


74 


M13A¥F2-3123 


-0.57±0.17 


6.()±1.5 


-238.3±1.0 


21 


M13A¥F3-548 


-0.49±0.16 


4.2±0.7 


-235.0±1.3 


26 


M13/WF3-1718 


+0.14±0.13 


2.0±0.3 


-243.7±0.7 


96 


M15/B203 


-0.14±0.22 


5.1±1.3 


-95.0±0.3 


32 


M15/B315 


+0.10±0.21 


1.8±0.2 


-103.2±0.9 


72 


M15/B334 


-1.38±0.86 


10.1±0.5 


-106.5±2.1 


11 


M15/B374 


+0.59±0.21 


4.0±0.9 


-107.0±1.1 


82 


M92/B176 


-1.26±0.75 


7.0±1.0 


-117.7±1.2 


21 


NGC288/B16 


+0.69±0.15 


2.9±0.9 


-45.1±0.9 


42 


NGC288/B22 


+0.52±0.22 


2.5±0.5 


-43.9±0.8 


47 


NGC288/B186 


+0.40±0.17 


3.3±0.8 


-43.9±0.5 


65 


NGC288/B302 


+0.77±0.17 


1.8±0.2 


-44.0±0.5 


75 



To process the spectra, Behr ( 2003ab has employed the package 
of routines developed by McCarthy ( 119901 ) for the FIGARO data 
analysis package i jShortridge 1993^ . A comprehensive description 
of the data acquisition and reduction procedure is presented by 
Behr ( I2003ab . 



3 SIMULATION PROCEDURE 

The line profile simulations were performed with the ZEEMAN2 
spectrum synthesis code (Landstreet ( 119881 ), Wade et al. ( I200U ). 
Khalack & Wade ('2006') have modified the ZEEMAN2 code to al- 
low for an automatic minimization of the model parameters using 
the downhill simplex method dPress et al. 1992t . To analyze the ver- 
tical abundance stratification we have built the dependence of the 
abundance derived from each analyzed profile relative to the opti- 
cal depth tsooo for the list of selected line profiles, assuming that 
the profile is formed mainly at line optical depth ^£=1. This depth 
corresponds to a continuum optical depth 1:5000 of a certain value 
that in turn corresponds to a given layer of the stellar atmosphere 
model. A comprehensive description of the procedure is given by 
Khalack et al. ( I2007t . This method was also used by Khalack et al. 
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( 12008 1 for the study of stratification in BHB stars and by Thiam et 
al. ( [20l0t for HgMn stars. 

The stellar atmosphere models used for our analysis were cal- 
culated with the Phoenix code (Hausch ildtetal. 1997} assuming 
LTE (Local Thermodynamic Equilibrium), solar metallicity except 
for the iron abundance and depleted helium abundances, and the 
surface gravities extracted from Behr ( |2003at which are listed in 
Table [T] The depleted helium abundance has also been taken into 
account during line profile simulations using ZEEMAN2. 

In a previous study, JKhalack et al. 2008l l showed that inclu- 
sion of non-zero microturbulence leads to the amplification of sig- 
natures of vertical abundance stratification for different chemical 
species. Therefore, in order to mitigate the influence of microtur- 
bulence on the estimation of iron vertical stratification, we have 
assumed zero microturbulent velocity in the simulations presented 
here. 

In this paper, we have selected a list of Fe II lines that are 
suitable for abundance and stratification analysis. The full list of 
the analyzed spectral lines is given in the electronic edition of 
MNRAS. Atomic data for Fe II lines are extracted from VALD-2 
(Kupka et al. jl999t , Ryabchikova et al. ('I999'l), Pickering et al. 
( 1200 It and Raassen & Uylingfl ( |1998| I) line database. 



4 IRON VERTICAL STRATIFICATION 

A sample of fourteen BHB stars from the globular clusters Ml 3, 
M15, M92 and NGC 288 were studied using the procedure briefly 
described in Sec.[3l These stars were chosen using the criteria de- 
scribed in Sec.|2]and some of their properties obtained by our anal- 
ysis are listed in Table |2l namely, their average iron abundance, 
V sin i and mean heliocentric radial velocity averaged over all lines 
analyzed (n represents the number of lines). The reported uncer- 
tainties are equal to the standard deviation calculated from the re- 
sults of individual line simulations for all lines considered. The re- 
sults shown in Table[2]for the iron abundance are consistent to those 
found by Behr ( .2003aJ (see Table[TJ except for the stars B334 and 
B176. This difference is partly due to the fact that Behr ( I2003at 
used mainly strong iron lines for his analysis while we also use 
weak lines which are formed deeper in the atmosphere and predict 
larger average iron abundance. A larger abundance is obtained be- 
cause these two stars show an enhancement of iron abundance in 
the deeper atmospheric layers. The radial velocities found here are 
also consistent for all the stars, but once again the stars B334 and 
B 176 show a small difference between our results and those of Behr 
( |2003a| >. These two stars will be discussed in more detail below. 

The effective temperatures of the BHB stars studied here range 
from 10,750K to 14,030K (see Tab.[B. This and the variety of the 
mean abundance values for iron (see Tab. |2} among the sample of 
BHB stars provides a unique opportunity to search for evidence of 
abundance vertical stratification of iron and to estimate its possible 
relation to other physical characteristics of these stars. 



4.1 Search for signatures of vertical stratification 

To verify if iron abundance is vertically stratified in the stellar at- 
mospheres of the selected BHB stars, we have analyzed individ- 
ually all Fell line profiles present in the spectra that give small 
errors for the estimated parameters (iron abundance, Vr and V sin;). 



Table 3. Data on vertical stratification of iron abundance. 



Cluster/Star 


Fell 


n 




slope 


l''g(^) 




M13AVF2-820 


0.30±0.02 


-4.59±0.32 


79 


M13AVF2-2692 


0.20±0.02 


-4.47±0.25 


74 


M13AVF2-3123 


0.05±0.07 


-5.07±0.17 


21 


M13AVF3-548 


0.04±0.04 


-4.99±0.16 


26 


M13AVF3-1718 


0.02±0.02 


-4.36±0.13 


96 


M15/B203 


0.15±0.04 


-4.64±0.22 


32 


M15/B315 


0.07±0.03 


-4.40±0.21 


72 


M15/B334 


0.72±0.08 


-5.89±0.87 


11 


M15/B374 


0.07±0.04 


-3.91±().21 


82 


M92/B176 


0.65±0.07 


-5.76±0.75 


21 


NGC288/B16 


0.01±0.03 


-3.81±0.15 


42 


NGC288/B22 


0.12±0.04 


-3.98±0.22 


47 


NGC288/B186 


0.13±0.02 


-4.10±0.17 


65 


NGC288/B302 


0.06±0.02 


-3.73±0.17 


75 



' I ftp://ftp.wins.uva.nl/pub/orthl 



assuming zero microturbulence. The measured abundance for each 
line was then associated to a line depth formation. The iron abun- 
dances obtained for different optical depths (see Figs. [Hi were 
fitted to a linear function (for the range of atmospheric depths 
—5.3 < logTsooo < —2.0) using the least-square algorithm to statis- 
tically evaluate the significance of observable trends. This range of 
logXgooo was chosen because spectra from our sample of stars all 
possess iron lines within this interval. The main result found here is 
that the iron abundance generally increases towards the deeper at- 
mosphere, although some of the increases found are not statistically 
significant. The slope of these gradients are presented in Table [3] 
where the individual columns give the increase of logAVe/Mo/ (in 
dex) calculated per dex of loglgoooi the average iron abundance 
value obtained taking into account all of the lines and the number 
of lines analyzed. 

Without further analysis, it is found that among our sample of 
fourteen stars, seven stars show clear signatures of vertical stratifi- 
cation of iron. Two of them (B22, B186) are located in the globular 
cluster NGC 288, two (WF2-820, WF2-2692) in M13, two (B203, 
B334) in M15 and one (B176) in M92 (see Fig.[I]and Tab.^Jl. One 
more BHB star, B302 (in NGC 288) also shows a slope that is sta- 
tistically significant (with respect to the error bar), but its value is 
too low to confidently state, due to all of the uncertainties inherent 
to line analysis, that iron stratification exists in this star 

However, for two of the seven BHB stars that show iron strat- 
ification (B334, B176), iron is strongly depleted (in comparison 
with the solar abundance of iron) and the number of analyzed iron 
lines is therefore respectively smaller than for most of the other 
stars studied here (see Tab. [3]l. However, the iron abundance at 
logXgooo = -2 is nearly solar and is therefore higher than the clus- 
ter abundance. The same two stars show very large slopes for the 
iron abundance with respect to optical depth (see bottom row of 
Fig.[T). Since our analysis is based on a relatively small number of 
lines (especially for B334), these results could be distorted by few 
cases of line misidentification. It is also surprising that such a large 
stratification slope is found where iron abundance is not enhanced 
in the outer atmosphere and helium abundance is large. If diffusion 
is at play, a larger iron abundance and a smaller helium abundance 
would be expected. These two stars also appear to have radial ve- 
locities, determined from the analysis of Fe II lines, slightly differ- 
ent (even though this difference is inside the error bars it is greater 
than for the other stars of our sample) than those obtained by Behr 
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Figure 1. Abundance estimates from the analysis of Fe II (filled circles) lines as a function of line (core) formation optical depth for BHB stars B203 in M15 
(middle top), WF2-82() (second row left) and WF2-2692 (second row right) in M13, B22 (third row left) and B186 in NGC 288 (third row right), B334 in M15 
(bottom left) and B176 in M92 (bottom right). The solid line approximates the data with a linear fit using the least-square algorithm. 
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( I2003ab from the combined analysis of lines of different chemical 
species. The average iron abundance found here for these two stars 
(see Tab. |2} are quite different form those found by Behr ( I2003ab 
(see Tab. [TJ, but this is due to the fact that we have also studied 
some weaker lines (see Sec. 14.1b . Therefore, all of these discrep- 
ancies cast some doubts about our results for these two stars and 
therefore more high-quality spectra is needed to confirm the pres- 
ence of iron stratification in their atmosphere. 

The vertical iron stratification for the other five stars: WF2- 
820 and WF2-2692 in M13, B22 and B186 in NGC 288 and B203 
in M13 is shown in Fig.[T] The lines used to produce these figures 
span several orders of magnitude for Xgooo- Also, since the increase 
of the abundance along this depth range is too large to be explained 
by the various uncertainties intervening in abundance determina- 
tion, it is safe to state that iron is vertically stratified in these stars. 



4.2 Dependence of stratification on effective temperature 

It is clear from the results shown in the previous section and those 
of Khalack et al. (2007 and 2008), that some BHB stars show sig- 
natures of vertical stratification of iron abundance, while others do 
not reveal such a property. All the studied stars (except B334, which 
has Teff ~10,750K) have Tgff higher than 1 1,000K and are therefore 
in or near the domain of Tgff where diffusion is expected to come 
into play. However, these stars are found in different globular clus- 
ters and are shown to have different mean iron abundances. The 
effective temperature is one of the most important parameter for 
a stellar atmosphere and has direct impact on the atomic diffusion 
in BHB stars (Hui-Bon-Hoa et al. 2000 and LeBlanc et al. 2009). 
It would therefore be very interesting to verify if a dependence of 
vertical iron stratification on T^.^ exists in BHB stars. 

To track possible dependence of iron vertical stratification on 
Teff we have combined the results obtained in this study with the 
previously published results by Khalack et al. ( 12008b for other BHB 
stars, whose iron abundance slope (see below) were recalculated for 
the same range of atmospheric depths — 5.3 < logXgooo < —2.0. 

The slope of the iron abundance as a function of optical depth 
for the various BHB stars was calculated and a correlation is found 
between this slope and r^ff. Fig. |2] shows the iron slopes with re- 
spect to Tgff, along with a linear fit for these data-points using 
the least-square algorithm (solid line). The stars B334 in M15 and 
B 176 in M92 were excluded from this fit because of the reasons dis- 
cussed in Subsec. |4ni They are however shown in Fig.|2]as empty 
circles. The inclusion of these two stars would not change the qual- 
itative conclusion that the slope of iron abundance decreases as a 
function of Teff for BHB stars with T^ff h 1 1,500 K. The BHB star 
WF3-1718 in M13 (shown as empty square in Fig.|2]l was also not 
taken into account during the aforementioned fit depicted by the 
solid line, because the wings of Balmer lines in its spectrum suggest 
that a higher Tgff is probable for this star. Nevertheless, inclusion of 
WF3-1718 still results in some visible correlation (see dashed line 
in Fig.|2]l- Adoption of higher effective temperature for WF3-1718 
would amplify this correlation. 

The detected correlation indicates that the strongest vertical 
stratification of iron abundance is expected for BHB stars with Teff 
near 11,500K while it decreases for higher values of Teff. It be- 
comes negligibly small for the BHB stars with Teff >z 14,000K. 
LeBlanc et al. (2010, in preparation) find that the theoretical models 
for BHB stars of LeBlanc et al. (2009) predict a similar tendency 
of decreasing iron vertical abundance with Tgff . 



5 DISCUSSION 

In this paper, we have found additional evidence that iron abun- 
dance is vertically stratified in the atmosphere of BHB stars. From 
the analysis of 14 BHB stars, selected as the most probable candi- 
dates for such detections due to their T^g y 11,500K and Vsini ^ 
10 km s^' , we have found clear signatures of iron stratification in 
five stars: B22, B186 in the globular cluster NGC 288; WF2-820, 
WF2-2692 in M13 and B203 in M15. Two other stars, B334 in M15 
and B176 in M92, show extremely strong increase of iron abun- 
dance as a function of optical depth in their atmosphere (see Tab. [3] 
and Fig.|2l(. However, this result may be influenced by the strongly 
depleted mean iron abundance (that results in a smaller number of 
analyzed line profiles). Also, because of other anomalies found for 
these two stars (see Sect. 14. lb . some doubts remain about observed 
stratification of iron in their atmosphere. 

For B334 and B176 the effective temperature is only slightly 
smaller than the limit (Tgg =11,500 K) above which BHB stars 
show the photometric jumps and gaps (Grundahl et al. (1999), Fer- 
raro et al. l |1998b ) and where diffusion is efficient. On the other 
hand, Khalack et al. J2008b have shown that stellar atmosphere 
models with higher Tf-g usually slightly decrease the slope of the 
element's abundance with respect to the logarithm of excitation po- 
tential as well as with respect to logXgooo. Therefore, adoption of 
higher effective temperatures for B334 and B176 will decrease the 
estimated slopes of iron abundance which increases towards the 
deeper atmosphere. 

The values of effective temperature for the sample of analyzed 
BHB stars were derived from Kurucz's ATLAS9 synthetic pho- 
tometry grids dBehr 2003ab . which use the stellar atmosphere mod- 
els with homogeneous distribution of element's abundance. Hui- 
Bon-Hoa et al. (2000) and LeBlanc et al. (2010, in preparation) 
have shown that the observable photometric jumps and gaps can 
be explained in terms of atomic diffusion mechanism that leads to 
stratification of element's abundance with atmospheric depths. The 
models with stratified abundance of chemical species result in dif- 
ferent photometric colors as compared to chemically homogenous 
models. This fact may introduce some errors in the evaluation of 
Teff by photometric means when using chemically homogeneous 
atmospheric models. The uncertainty on Tgff determination should 
be studied by using the new atmospheric models of LeBlanc et al. 
(2009) that take into account vertical stratification of the elements. 
This is outside the scope of the present paper. 

Our analysis of the dependency of the available data for verti- 
cal iron stratification in BHB stars with respect to the effective tem- 
perature results in the detection of a trend (see Fig. |2}. The same 
results were obtained from the modeling of iron vertical stratifica- 
tion employing the self-consistent stellar atmosphere models with 
vertically stratified abundance of chemical species (LeBlanc et al. 
2010, in preparation). The results presented here show that for BHB 
stars with Tgff around 1 1,500K have the strongest vertical iron strat- 
ification gradients. This gradient diminishes as Teff increases up to 
Teff ~ 14,000K, where no significant iron stratification exists. Nev- 
ertheless, to statistically improve the registered trend of iron strat- 
ification more data for BHB stars with appropriate values of Teff 
and Vsini, especially for Tgff near 11,500K need to be analyzed. 
Such additional data will be useful to constrain theoretical model 
atmospheres that include stratification of elements, similar to those 
of LeBlanc et al. (2009). 
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Figure 2. Plot of the slope of iron vertical stratification versus effective 
temperature of analyzed BHB stars (see text for details). The dashed line 
is a linear fit that includes all stars except for the two stars represented by 
empty circles (B176 and B334). The linear fit given by the solid line also 
excludes the star WF3-1718 represented by the empty square. 
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